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Abstract 

'NT 

We report on simulations of DT simplicial gravity for manifolds with the topology of the 
' 4-disk. We find evidence for four phases in a two-dimensional parameter space. In two of 

I these the boundary plays no dynamical role and the geometries are equivalent to those observed 

I earlier for the sphere 5^. In another phase the boundary is maximal and the quantum geometry 

' degenerates to a one dimensional branched polymer. In contrast wc provide evidence that the 



X 



fourth phase is effectively three-dimensional. We find discontinuous phase transitions at all the 
phase boundaries. 



Introduction 

Dynamical triangulation (DT) models arise from simplicial discretizations of continuous Rieman- 
nian manifolds. A manifold is approximated by glueing together a set of equilateral simplices with 
fixed edge lengths. This glueing ensures that each face is shared by exactly two distinct simplices - 
the resultant simplicial lattice is often called a triangulation. In the context of Euclidean quantum 
gravity it is natural to consider a weighted sum of all possible triangulations as a candidate for 
a regularized path integral over metrics. Physically distinct metrics correspond to inequivalent 
simplicial triangulations. This prescription has been shown to be very successful in two-dimensions 
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in which analytic methods have been complemented by simulation studies, (see, for example, Q). 
In dimensions above two there are no known analytic techniques for handling the sum over tri- 
angulations and we must rely on numerical simulation. All of the latter studies have focused on 
elucidating the phase structure of compact manifolds, principally the sphere S"^. In this paper we 
investigate the phase structure of four manifolds with the topology of a 4-disk - that is a four sphere 
5"^ equipped with a single boundary with topology S^. As we will demonstrate the dimension of 
the parameter space of this model is larger than the corresponding compact models. The simplest 
compact models exhibit only discontinuous phase transitions precluding a continuum limit. One 
motivation for the current work was to see whether the richer parameter space of the non-compact 
models contains any continuous phase transitions. A natural lattice action Sb can be derived from 
the continuum action by straightforward techniques Q. It contains both the usual Regge curvature 
piece familiar from compact triangulations together with a boundary term. The boundary term 
arises from discretization of the extrinsic curvature of the boundary embedded in the bulk. In 
four-dimensions the curvature is localized on triangles. If Tm denotes the set of triangles in the 
bulk of the 4-triangulation (excluding the boundary) and Tqm those in the boundary the action 
can be written 

Seh = /«2 ~ ^"^h) + X! ~ ""''^) + ^4^'4 (1) 

\heTM heTsM J 

The quantity a = arccos (1/4), is the 4- volume and Uf^ is the number of simplices sharing the 
triangle (hinge) h. The curvature part of the action can be rewritten in terms of the number of 
vertices Nq, the boundary volume A'^^ and the number of boundary vertices Aq. With this in mind 
we shall consider the general simplicial action 

Sb = -KoNo + K4A4 + KbNl + (2) 

This form of the action contains both cosmological constant terms and curvature terms for the 
bulk and boundary. Notice that the term involving is absent in three dimensions since it is 
related to the curvature of a boundary two-sphere which is a topological invariant. In this work we 
have always set = and K4 is used to tune the volume of the 4-disk. We are thus left with a 
two-dimensional phase space parameterized by kq and Kb conjugate to the number of vertices and 
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the number of boundary tetrahedra. 

The partition function for the system is then 

Z = Y: e~'^ (3) 

T 

where the sum is over triangulations, T. 

Simulation 

Our simulation algorithm is an extension of the algorithm for compact manifolds in arbitrary 
dimension described by Catterall P] and is described in [11|. To simulate a triangulation with 



a boundary we actually simulate a compact triangulation with the topology of the sphere but 
consider one marked vertex to lie outside the bulk triangulation. This vertex may never be removed 
during the course of the simulation and the surface of the ball created by its neighbour simplices 
constitutes the triangulated boundary of the 4-disk. The form of the action (equation |2|) can then 
be derived by writing down an expression for the integrated curvature of the full sphere expressed as 
contributions from the 4-disk and the simplices around the marked vertex. The latter contributions 
can be evaluated explicitly in terms of bulk and boundary simplex numbers and yield this simple 
form for the action. To perform a simulation we need a set of local moves which are ergodic on the 
space of triangulations of the 4-disk. Fortunately we have such a set of moves - the usual moves on 
the sphere. 

In four-dimensions there are just 5 types of move: vertex insertion, vertex deletion, exchange of a 
link with a tetrahedron (two moves: link to tetrahedron and tetrahedron to link), and exchange 
of one triangle for another triangle. Where these moves take place on sections of the triangulation 
involving the marked vertex we take care to count changes in the numbers of simplices inside and 
outside the boundary so that we can calculate the change in the action, but otherwise the moves 
are the same as for the bulk. 



The code is described in more detail in |11] where it was used for the simulation of three-dimensional 
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dynamical triangulations with a boundary. The code was written for arbitrary dimension and earher 
checked against other workers' results in two-dimensions O. 



We have used the Metropolis Monte Carlo [10| scheme with usual update rule: 

^(accept move) = min{e^^^'' , 1} (4) 
and in this way we explore the space of triangulations with the action Sf,. 

Measurements 

Here we define some of the measurements used to characterize the configurations obtained from our 
simulations. During the simulations we store configurations at some interval which is sufficient to 
ensure that the configurations are independent. We check this by estimating the auto-correlation 
time (r) of each measurement when calculating expectation values. Uncorrelated data would give 
r = 0.5, we quote r when it significantly exceeds 0.5. 

We use two geodesic measures, davg and dhdy- The average geodesic distance between simplices, 
davg, is measured by counting the smallest number of steps between adjacent simplices required 
to get from one randomly selected simplex to another randomly selected simplex and taking the 
mean over a sample of such measurements. The mean geodesic distance to the boundary, dbdy, 
is measured by counting the smallest number of steps between adjacent simplices required to get 
from one randomly selected simplex to the boundary and taking the mean over a sample of such 
measurements. The last step from a boundary simplex to the boundary is counted as 0.5. 

To discuss the singular vertex structure observed in some phases we will consider a coordination 
measure, q, defined as the number of simplices sharing a given vertex. Similarly we define (qmax) 
as the expectation value of the coordination of the most highly coordinated vertex; (qnextMax) as 
the expectation value of the coordination of the next most highly coordinated vertex; {qavg) as the 
expectation value of the mean coordination of all vertices in the triangulation; and (qbdy) as the 
expectation value of the coordination of the boundary by which we mean the average number of 
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unique 4-simplices being shared by a typical boundary vertex. 

Phase diagram 

We performed a set of simulations in four dimensions with action of equation |2|. In all runs K4 was 
used to tune the nominal system volume, for each given kq and k^. To map the phase diagram 
we used A'4 = 1000; while to characterize the phases we used A^4 = 8000. In order to check the 
orders of transitions we have also used finite size scaling using additional simulations at A^4 = 2000 
and A^4 = 4000. 

Series of runs varying either hq or Kb were made and the vertex susceptibility used to search for 
phase transitions. We define the vertex susceptibility, x> to be normalized with respect to the 
number of 4-simplices: 

X = ^((iVo') - (Nof) (5) 

The points shown in figure Q are taken from the positions of peaks in the vertex and boundary 
susceptibilities. 

In figure ^ there are four phases which we characterize as: crumpled, minimal boundary (CMB); 
branched-polymer, minimal boundary (BPMB); boundary dominated (BD); and intermediate bound- 
ary (IB). 

In CMB and BPMB phases the boundary is simply 5 tetrahedra (3-simplices) connected to form 
a hyper-tetrahedral hole. The system is essentially a sphere with one marked 4-simplex — the 
hyper-tetrahedral hole. 

CMB phase 

Here we show typical data for ^^^^4 = 8000 at kq = —4 and k^, = 4 using 137 samples (a sample corre- 
sponds to 100000 attempted updates). The boundary size is just 5 - a minimal hole and we see one 
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Figure 1: Boundary size ((iV|)) and number of vertices {{Nq)) for 4-dimensional dynamical trian- 
gulation with a boundary. Nominal simulation volume, A'4 = 1000. 

'singular' vertex with (qmax) = 2642(6) (the next most coordinated vertex has {qnextMax) = 531(5), 
and {qa^g) = 250.7(5)) 

The presence of singular vertices is similar to the behaviour seen in compact triangulations in 4- 
dimensions [Q, ^. In the compact case, there are two singular vertices which are equally coordinated. 
In simulations with a boundary we find one singular vertex and the minimal-boundary appears to 
assume the role of the other singular vertex (the number of simplices shared by boundary vertices 
is 1871(82)). 

We find {dMy) = 9.18(4) (some correlation, r = 1.7) and (davg) = 10.66(1). The approximate 
equality of these two distances tells us that in this phase the presence of the boundary plays no 
crucial role (except for changing the number of singular vertices). Furthermore, we have observed 
that the typical manifold is very crumpled, having only a logarithmic growth of its mean size. Such 
a behavior may be characterized by a large effective dimension and is reminiscent of the usual 
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Figure 2: Phase diagram for 4-dimensional dynamical triangulation with a boundary. Nominal 
simulation volume, = 1000. Error bars are from estimation of the positions of the susceptibility 
peaks. Where error bars cannot be seen they are smaller than the symbols; the lines are guides to 
the eye. 

crumpled phase seen in simulations of the sphere. 



BPMB phase 



Again we characterize the phase using data for = 8000 (116 samples have been acquired) with 
couplings kq = 5, Kb = 6. We find {dbdy) = 43.3(10) and {davg) = 45.6(5). Again, the distance 
to the boundary is comparable to the distance to any randomly selected simplex indicating that 
the former has no distinguished role in the triangulation. Furthermore, large mean-geodesics are 
consistent with paths constrained to follow long branches. In this phase the vertex coordinations 
show no sign of singular structure: {qmax) = 387(7) and the next most highly coordinated vertex 
is about 50 lower, (qavg) = 20.90(1). Thus, the physics of this phase is again rather like the 
corresponding situation on the sphere, the typical geometries are one-dimensional polymers. 
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BD phase 



In this phase the boundary is essentially maximal, the maximal value being obtained from an 
arbitrarily branched chain of simplices such that N^maximal = 3 x Ar4 + 2. This is easy to understand 
in terms of a single chain of simplices - the 2 end 4-simplices have 4 boundary 3-simplices, and 
the (A^4 — 2) middle 4-simplices each have 3 boundary 3-simplices. If one then considers adding a 
branch, one boundary simplex is lost at the branch point but one boundary simplex is gained at 
the end simplex, so there is no net change in N^. 

To support this we list data from a run at A'4 = 8000 and couplings kq = —4, ki, = —10, in which 
we have collected 167 samples. We find that the boundary size fluctuates only very slightly about 
a mean of 7V| 24002 = 3 * 8000 + 2 and is very branched, possessing (A^end) = 2534(2) end 
points. The maximum vertex coordination {qmax) = 103(2) is close to that expected for a flat 
lattice. Another indicator that the geometry corresponds to narrow tubes is seen when we examine 
the geodesic distances: {dbdy) = 0.5038(2) and (davg) = 125.8(2) The mean boundary distance 
tells us that every simplex is a boundary simplex. Thus the quantum geometry is effectively one- 
dimensional. We find {qavg) = 4.9975000(1) which we take as further indication that there are no 
sections of bulk in this phase. The argument for this is as follows: in a minimal-width chain with 
no branches we have two ends which have a few vertices shared by less than 5-simplices — the 
coordination of vertices in the bulk of the chain. At each end there is 1 vertex with coordination 
1, 1 with coordination 2, 1 with coordination 3, and 1 with coordination 4. This adds up to a 
coordination deficit of 10 for the 4 vertices at an end, 20 for both ends. In a minimal-width chain 
consisting of N4 simplices we expect Nq = N4 + A and so with N4 = 8000 we expect, and see. 
No = 8004. This allows us to calculate the expected (qavg) = (iVo x 5 - 20)/iVo, putting iVo = 8004 
wc get {qavg) = 4.997501. Adding branches to the chain results in more highly coordinated vertices 
at the branch point that exactly cancel the deficit from the additional end. 
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IB phase 



Here we show data for A'4 = 8000 (86 samples) at couplings kq = —4, Kb = 0. We find that, once 
again, the boundary size scales linearly with the volume (see figure]^ yielding ~ x 1.033. We 
also find one large singular vertex, with {qmax) = 5721(9) {{qnextMax) = 233(2), (qavg) = 133(1)). 
The measurements of geodesic distances show that the boundary plays a preferred role in the 
triangulation (dbdy) = 0.834(2) while (davg) = 11.36(1). These numbers are incompatible with 
an extremal branched polymer shape. However, these measurements by themselves would not be 
inconsistent with a 'fat-branch' model. If we add to this the presence of a single singular vertex and 
the observation that about 40% (3495 (9) out of 8000) of the simplices have one boundary face we 
tentatively conclude that the typical geometry in this phase is three-dimensional — the boundary 
of the system coinciding with the boundary of the 4-ball surrounding the singular vertex. 

100000 r ■ ' ' ' ■ ■ I 




1000^ . : : : • • ' ' ' ' 

1000 10000 100000 

Simulation volume, N4 



Figure 3: Plot showing scaling of boundary size with simulation volume for the intermediate 
boundary (IB) phase, = —4, Kh = 0. The error bars are much smaller than the symbols. 
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Phase transitions 



The phase transition between the CMB and BPBM phases has been studied in the simulations of 
compact systems in three and four-dimensions and found to be discontinuous (first-order) in both 
cases (3d 4d ^). Our simulations have minimal boundary in both phases and so we expect 
the same behaviour. This was verified in three dimensions [pi] ]. 

We have also investigated the transition between CMB and IB phases and find good evidence here 
also for a discontinuous phase transition (figure top; N4 = 1000, kq = —4, ki, = 0.9). Similarly, 
the time series close to the BPMB — BD phase boundary (figure ^ middle; A^4 = 1000, kq = 4, 
Kfe = 1.64) also shows signs of bistability indicative of a discontinuous transition. 

At = 1000, 2000 and 4000 we found no such signals near the BD — IB boundary. However, we 
found linear scaling of the height of the peak in the boundary-size susceptibility which indicates 
a first-order transition. To confirm this we extended our simulations to N4 = 8000 and found 
bistability in the time series (figure bottom; A'4 = 8000, kq = —4, Kb = —1.862). 

Concluding remarks 

We have simulated four-dimensional simplicial gravity on manifolds with the topology of a 4- 
disk. Our action contains both bulk curvature terms and a boundary cosmological constant term. 
We have identified four phases in the model within the range of couplings —6 < hq < 6 and 
—4 < Kf) < 8. The observed phases include the crumpled and branched-polymer phases seen in 
triangulations of compact manifolds, and the boundary dominated phase seen in three-dimensions. 
The latter consists of a maximally branching tree. We have also identified a fourth phase which 
resembles a singular 4-ball in which essentially all simplices share a common bulk vertex and the 
physics is dominated by the three-dimensional boundary. 

All of the boundaries between these phases appear to be associated with discontinuous phase 
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Figure 4: Simulation time series showing bistability in the boundary size (-/Vf) at the CMB — IB 
phase boundary (upper; N4 = 1000, kq = —4, = 0.9); the BPMB — BD phase boundary (middle; 
N4 = 1000, Ko = 4, K5 = 1.64); and the BD— IB phase boundary (lower; iV4 = 8000, = -4, 
Kb = —1.862). We take this as indication of the discontinuous nature of these transitions. 

transitions. This is the same situation as for four-dimensional simplicial gravity on manifolds with 
topology and means that we cannot take a continuum limit in the vicinity of any of the observed 
phase transitions. Notice that all four phase boundaries appear to meet (within our errors) at one 
unique point. The simplest explanation for this feature is to assume that our phase diagram 
contains only two independent transition lines which would then generically intersect at a single 
point. Glancing at the nature of the phases one would associate one line with a boundary dividing 
lattices with singular vertices from those without. The second line would then correspond to a 
dividing line between geometries with extended boundary and those with minimal boundary. 

Notice also, that in all the phases we have identified, the number of vertices on the boundary is 
strongly correlated with the boundary volume. This need not be the case and it would be interesting 
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to vary the coupling away from zero to see whether a phase could be found which exhibited a 
classical scaling of boundary size with 4- volume. 

Acknowledgments 

Simon Catterall was supported in part by DOE grant DE-FG02-85ER40237. Simeon Warner was 
supported in part by the DOE under contract W-7405-ENG-36. 

References 

[1] E Adi, M Hasenbusch, M Marcu, E Pazy, K Pinn and S Solomon. Monte carlo simulation of 
2-d quantum gravity as open dynamically triangulated random surfaces. Physiscs Letters B, 
320 227-233 (1994). arXiv:hep-lat/9310016. 

[2] J Ambjorn. Quantization of geometry. (1994). arXiv:hep-th/9411179. 

[3] J Ambjorn, D V Boulatov, A Krzywicki and S Varsted. The vacuum in three-dimensaional 
simplicial quantum gravity. Physiscs Letters B, 276 432-436 (1992). 

[4] P Bialas, Z Burda, A Krzywicki and B Petersson. Focusing on the fixed point of 4d simplicial 
gravity. Nuclear Physics B, 472 293-308 (1996). arXiv:hep-lat/9601024. 

[5] S Catterall. Simulations of dynamically triangulated gravity - an algoritm for arbitrary di- 
mension. Computer Physics Communications, 87 409-415 (1995). arXiv:hep-lat/9405026. 

[6] S Catterall, R Renken and J Kogut. Singular structure in 4d simplicial gravity. PhysLettB, 
416 274-280 (1998). arXiv:hep-lat/9709007. 

[7] S Catterall, G Thorleifsson, J Kogut and R Renken. Singular vertices and the triangulation 
space of the d-sphere. Nuclear Physics B, 468 263-276 (1996). arXiv:hep-lat/9512012. 



12 



[8] B de Bakker. Further evidence that the transition of 4d dynamical triangulation is Vt order. 
Physiscs Letters B, 389 238-242 (1996). arXiv:hep-lat/9603024. 

[9] J B Hartle and R Sorkin. Boundary terms in the action for the regge calculus. General 
Relativity and Gravity, 13 541-549 (1981). 

[10] Nicholas Metropolis, Arianna W Rosenbluth, Marshall N Rosenbluth, Augusta H Teller and 
Edward Teller. Equation of state calculations by fast computing machines. The Journal of 
Chemical Physics, 21 (1953). 

[11] Simeon Warner, Simon Catterall and Ray Renken. Phase diagram of three-dimensional 
dynamical triangulations with a boundary. Physiscs Letters B, 442 266-272 (1998). 
arXiv:hep-lat/9808006. 



13 



